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Abstract 

We report experimental implementation of discrete Fourier trans- 
formation(DFT) on a nuclear magnetic resonance (NMR) quantum 
computer. Experimental results agree with theoretical results. Us- 
ing the pulse sequences we introduced, DFT can be realized on any 
L-bit quantum number in principle. 
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1 INTRODUCTION 

With the discovery of an apparent separation between the classical and 
quantum classifications of computational complexity[|l]],and of fault-tolerant 
schemes for quantum computation [Q, quantum information theory plays an 
important role in computer science. In 1994, Shor|jl]] introduced the quantum 
factoring algorithm, which achieves an exponential speed-up relative to clas- 
sical algorithms. Shor's crucial insight was that the discrete Fourier transfor- 
mation(DFT) can be evaluated in polynomial time on a quantum computer. 
In order to realize Shor's algorithm, performing DFT becomes most impor- 
tant. In this paper, we introduce a version to realize DFT by using NMR 



spectrometer and simulator and give our experiment and simulation results 
realizing two qubit DFT. 

The discrete Fourier transformation modulog(DFTg) is a unitary trans- 
formation in g-dimension(g = 2^). It is defined relative to chosen basis 
|0),---,|g-l)byi: 

DFTg :| a) — ^ ^ exp{2mac/q) \ c) (1) 

v'' c=0 

Coppersmith [|] suggested the efficient algorithm for DFT on application 
of quantum mechanical operators and gave an improvement in which evalu- 
ation of L bit Fourier transformation is accomplished by composing L one- 
qubit operations and — 1) two-qubit operations. So the implementation 
of DFT is based on realizing one-qubit operator Aj which acts on the state 
of the j-qubit and two-qubit operator Bj^ which acts on the states of qubit 
j and k. The explicit form of the operator Aj is[0: 

A, = 2-1/2(1 0^.^(0, I + I 0,)(1, I + I 1,)(0, I - I I) 

(J = ) 

The matrix representation of is: ^ ^ | ^ 

The explicit form of the operator Bjk is[^: 

B, k =1 0,0,)(0,Ofc I + I 0,1,) (0,U I + I ljOfc)(l,Ofc 

(3) 



(2) 



The matrix representation of Bj^ is: 



/ 1 



V 



Transformation of Bjk affects only state 

We use a network to illustrate how to use operators Aj and Bjk to compose 
5-qubit transformation in Figure. 1. The operators corresponding to the 
network are: 

(Aa) (-834^43) (524-823^2) (-Bi4-Bi3-Bi2^i) (-Bo4-Bo3-So2-Boi^o) 



2 IMPLEMENTING THE DFT IN NMR 

In this paper we realize DFT by using NMR. The prime thing we should 
do is to perform Aj and Bjk by using radiofrequency pulses and spin-spin 
interaction. Using NMR to perform quantum computer, we should choose 
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AX couple system. The Hamiltonian for this system can be approximated 
as: 

H = 2wabJzaJzB + '^aJzA + ^bJzB + Henv (4) 

We denote the rotation v^x'^^ = exp(i(/?/2,.) for a Lp rotation about the x-axis 

and ify^'^ = exp{iiply) for a ip rotation about the l/-axis, j or k means which 
nucleus is operated on, and r is the time within which the system undergoes 
the unitary transformation exp{2TTiJzAJzBt) in the doubly rotating frame. 

To perform Aj is actually to perform a Walsh-Hadamard transformation, 
which rotates each quantum qubit from | 0) to (| 0)+ | l))/-\/2. the pulse 
sequences is|[: (f)x - (f )x - (-f )y 

To perform Bj^, we use radiofrequency pulses and spin-spin interaction 
like: 

(-1 yy(-f - i-^yxj-vYx - m&y - i - (^)x - i - i^Yxi^rx 

Using the pulse sequences and spin-spin interaction we introduced above, 
we can implement DFT to any quantum number with L-qubit in principle. 



3 EXPERIMENTAL AND SIMULATION RE- 
SULTS 

In order to demonstrate the results described above, we have constructed an 
NMR quantum computer capable of implementing the DFT. In this section, 
we give NMR experimental and NMR simulation results after doing DFT on 
state I Ojlfc) by using the two radiofrequency pulse sequences corresponding 
to the operator Aj and Bjk respectively. Data of the NMR experiment were 
taken at room temperature with a Bruker ARX-500 spectrometer and the 
signals were obtained with single-shot measurement. 

In the experiment we chose H2PO3 as our sample, labeled ^^P as j- 
qubit and as fc-qubit. The observed J-coupling between and ^^P was 
64:7.4:51Hz. First, we produced effective pure state | 00)by using "temporal 
averaging' [H . The pulses we used were|0]: 

E: NONE 

r>i . (TT\P 1 T (■rT\p (■rT\H 1 (Tt\H 
ri. y^lX 2"' \2)Y\2)X 2J \2)Y 
TDH). (Tr\H 1 7 (n\H (iT\P 1 ('^\P 
\2)X 2"' \2)Y\2)X 2J ^ifY 

Since we wanted to perform DFT on the state| 01), so after the prepara- 
tion of the state | 00), we should operate the pulse sequence: (f )y — jj — (f )x 
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on to obtain the state| 01). 

Second, we performed discrete Fourier transformation on 
the pulses and spin-spin interaction are as follows: 



01). We list 



X 



-i-m- 



'n\P {n\H 



.JJXKJJX ' 



The first two pulses and the last two pulses correspond to the operator Aj 
and Ak respectively, the rest pulses correspond to Bj^. After these operations, 
the initial state| 01) was transformed into state: 

I (I 00)- I 01) +2 I 10) -i I 11). (5) 

Compared with the definition of DFT (1) which means that after doing 
DFT operation, the state | 01) should become 

I 01) — >l (I 00) + I 01) + e'"" \ 10) + e''^ \ 11)) 



(I 00) + z I 01)- 



and the density matrix is: 



(7) 



10) -i I 11)) (6) 
/I -i -1 i \ 
i 1 -i -1 
-1 z 1 -i 
\i -1 i 1 J 
The third step: we should reverse the qubit. That is, for example, for 
the case of three qubits, | ijk) — kji). The third step is composed of three 
C-NOT operations. The network is shown in Figure. 2. The radiofrequency 
pulses corresponding to the three C-NOT operators are: 



3. 
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Actually, the quantum-mechanical operation "reverse the qubit" is not 
applied. Instead of this, one measures the state after DFT and reads the 
result of the measurement in the opposite order ||^. 

After these three steps, we had finished the operation of DFT on the 
initial state | 01) . In order to illustrate the result we got, we need to obtain 
all the elements in the two-spin density matrix of the ultimate state by using 
state tomography^. Before we did the experiment, we had simulated the 
experiment by using NMR simulator. The results are shown in Figure 3a , 
it can be regarded as a theoretical result. The matrix representation of the 
/I -i -I i \ 

which is the same as we expected. 



result is: 



-1 



V 



1 

-i -I i 1 

An advantage of doing NMR simulator is that it can represent the physics 
problem in spite of considering experimental errors (chemical shift, the affec- 
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tion of environmental and so on). So using NMR simulator can be considered 
as a quantum computer to simulate some quantum mechanical problem. 

The experimental results are shown in Figure 36, we measured and got 
real and imaginary components of the integral area from each peak, then 
we used all the integral areas we'd got to reconstruct that density matrix. 
Compared with the theoretical results, the experimental results agree with 
those deduced from theory. The relative errors are 24.8% and 4.7%. They 
primarily due to the imperfect calibration of the rotation, and are also caused 
by inhomogeneity of the magnetic field or magnetization and least-squares 
fitting used in the tomography procedure. Using other NMR techniques, 
such as phase cycling, the errors can be minimized. 



4 CONCLUSION 

We have demonstrated that the discrete Fourier transition can be imple- 
mented by means of NMR spectrometer and simulator and the experiment 
agrees with the theory well. Using the method we introduced, we can perform 
DFT on any quantum number with L-qubit. Performing DFT successfully 
gives a probability to realize Shor's algorithm by using NMR. Fourier analy- 
sis is a versatile tools in the laboratory |T^, so we might expect that the fast 
DFT should be an important application to physics. 
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Figure 1: Quantum network for the quantum discrete Fourier transform 



Performing on a five-qubit register, gate B^^ performing the conditional phase 



shift by 2 Ji 



B 



Figure 2: The network of performing operation "reverse the qubit" 
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Figure 3. Theoretical(simulator) and experimental(spectrometer) deviation 
density matrices after operated DFT on state |01>. The unit of magnitude is 
arbitrary and the unit of phase is n . Relative errors were abtained 

according to : | | P theory" P expt I I / I I P theory I 



